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The dynamical evolution of the scale factor of FRW cosmological model is presented, when the 
equation of state of the material content assume the form p = 7/3, 7 = constant, including the 
cosmological term. We use the WKB approximation and the relation with the Einstein-Hamilton- 
Jacobi equation to obtain the exact solutions. 

PACS numbers: 04.20.-q, 04.20.Jb, 98.80.-k 



I. INTRODUCTION 

The behaviour of the cosmological scale factor A(t) 
in solutions of Einstein's field equations with the 
Friedmann-Robertson- Walker line element has been the 
subject of numerous studies, where the presentations 
tend to focus on models in which p = and there is 
no cosmological constant (A = 0). Some treatments in- 
clude the cosmological constant Q, Q> El S Q an d the 
pressure p is given in terms of density p by an equation 
of state p = pip) and 0, for particular values in the 
7 parameter |d la Isl ■ 

The standart model of cosmology is based on Einstein's 
General Relativity theory, which can be derived from the 
geometric Einstein-Hilbert Lagrangian 
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gco 



16ttG 



-gR, 



(i) 



where R is the Ricci scalar, G the Newton constant, and 
g = IgjLwl the determinant of the metric tensor. By per- 
forming the metric variation of this equation, one obtains 
the well known Einstein's field equations 



Raw — 2^£m^ = —8ttGT^ u , 



(2) 



where T^„ is the energy-momentum stress tensor, as- 
sociated with a matter lagrangian, which is the source 
of gravitation, assigning the corresponding equation of 



state, which varies during different epochs of the history 
of the universe. 

Introducing a symmetry through the metric tensor, in 
cosmology one assumes a simple one according to the cos- 
mological principle that states the universe is both homo- 
geneous and isotropic. This homogeneous and isotropic 
space-time symmetry was originally studied by Fried- 
mann, Robertson, and Walker (FRW). The symmmetry 
is encoded in the special form of following line clement 



ds 2 



-N 2 (t)dt 2 + A 2 (t) 



dr 2 



1 - KV< 



r 2 df7 z 



(3) 



where A(t) is the scale factor, N(t) the lapse function, 
k is the constant curvature, taking the values 0, +1,— 1 
(flat, closed and open space, respectively). 

The FRW solutions to the Einstein field equation 
represent a cornerstone in the development of modern 
cosmology, since with them it is posible to undertand 
the expansion of universe. 

Recently, Faraoni 10] introduced one procedure based 
in the Riccati differential equation, obtained by the 
combination of the Einstein field equation, resulting in 
the same solutions obtained by the standard procedure 
D IE B S IE E3, without the cosmological term. This 
alternative approach is more direct than the standard 
one, which was used in the factorization procedure in 
the supersymmetric level ^E 13 f° r to obtain the both 
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iso-spectral potential and function in particular one di- 
mension systems. 

The set of differential equation for the FRW cosmolog- 
ical model, including the cosmological term, become 
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(4) 
(5) 



the overdot means d/dt. 

In the literature one can find the well-known classical 
behaviour to the scale factor for k = \v\ 



A = 



6ttGm 7 ( 7 + \y 



(t - f ) 3h+1) • 



(6) 



Taking different values for the constant 7 we have the 
following subcases 



A = 



f irGMi 

o 3 



t* 



[6vrGM ] 3 ^ 
[24ttGMi]M 

However, for the case 7 = 
exponential 



(7) 



for 7 = i radiation 

for 7 = dust 
for 7 = 1 stiff fluid 

— 1 the solution becomes to 



A = A e 



Ht 



with H = 2\ -7rGM 



(8) 



here, we consider the sign (+) in the exponential func- 
tion, because we consider the inflationary behaviour. 
These solutions will be compared with the solutions 
found by our method. 

The main purpose of this work is the introduction of 
the WKB-like procedure for to calculate the function 
A(t), including the function $>(A), function that play 
a important role in the supersymmetric fashionjbll 
called the superpotential function, into the Hamiltonian 
formalism for to solve the Einsten-Hamilton- Jacobi equa- 
tion. Also, we include the cosmological term in the for- 
malism. 

The remainder of the paper is organized as follow. 
The procedure that include the Einstein-Hamilton- Jacobi 
equation and the master equation is described in Sec. II. 
In Sec. Ill we present the exact solutions for the master 
equation found for this model, including the correspond- 
ing analysis of them. Finally, the Sec. IV is devoted to 
comments. 



II. EINSTEIN-HAMILTON- JACOBI EQUATION: 
THE WKB-LIKE METHOD 



We will use the total Lagrangian for a homogeneous 
and isotropic universe (FRW cosmological model), and 
perfect-fluid like ordinary matter with pressure p and 
energy density p, and barotropic stat e eq uation p = 7/9, 



including the cosmological term A 



itate eqv. 



fiA /d A \ 

L= — f— J — 6kNA— 2N A A 3 + 167rGNM 7 A~ 37 . (9) 

We define the canonical momentum conjugate to the gen- 
eralized coordinate A (scale factor) as 

<9L f 2 A dA 



n A 



(10) 



dA N dt ' 
The canonical hamiltonian function has the following 
form 



L = n A A - NH = U A A 



-N 



TT 2 

24A 



6kA + 2AA* 



16ttGM 7 A" 37 



(11) 



where 



ft = 777T [n| + I44kA 2 + 48AA 4 - 384ttGM 7 A- 37+1 1 . 
24A L A 

(12) 

Performing the variation of (| IL 1 1> with respect to N, 
3L/9N = 0, implies the well-known result 7i = 0. 

AT this point we can do two things: i) The quan- 
tization procedure, imposing the quantization condition 
on H — > H, where Ji is an operator, and applying this 
hamiltonian operator to the wave function 'J, we obtain 
the Wheeler-DeWitt (WDW) equation in the minisuper- 
space 



= 0, 



(13) 



and ii) WKB like method, if one perform the transfor- 
mation 



TT d$ 



(14) 



in (|12|) . becomes the Einstein-Hamilton- Jacobi equation, 
when $ is the superpotential function that is related to 
the physical potential under consideration. 
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We shall use the part ii) as an alternative method for 
obtain the classical solutions to the FRW cosmological 
model. 

Introducing the ansatz p4|l into the Eq. (|12f> we get 



^ ) + 144kA 2 + 48AA 4 - 384 7 rGM 7 A- 3 ' , ' +1 
dA ' 



d<E> 8 A 

— = ±12AW-7tGM^A-( 3 '/+ 1 ) - — A 2 - k. 
dA V 3 7 3 



= 



(15) 



Relating the equations I|10ll4fl and l|15[) . we obtain the 
classical evolution for the scale factor in term of the "cos- 
mic time" r defined by dr = N(t)dt, through the follow- 
ing master equation 



dr : 



dA 



| 7 rGM 7 A-(^+ 1 ) - f A 2 - k 



(16) 



that correspond to Eq. (0 in the gauge N=l. 

This equation is not easy to solve in general way for all 
values in the 7 parameter. However, we can solve this one 
for particular values in two sectors in the 7 parameter: 

1. 7 < 0, it it say (-1/3, -2/3,-1) and A ^ 0. This is 
the phenomenon commonly known as inflation-like. 

2. 7 = 1/3, A^O, any k. 

3. 7 > 0, A = 0. 

In the following, we describe its behaviour for the scale 
factor. 

III. SOLUTION TO THE MASTER EQUATION 

Here, we obtain analytic solutions for the scale factor, 
via the master equation rewritten in terms of a "confor- 
mal time" coordinate T. In some cases, will be necessary 
to drop the cosmological term, with the end of obtaining 
the corresponding exact solution. 

A. 7 < 0, inflation-like phenomenon 

Considering some negative values for the 7 parameter, 
namely, 7 = —1, —1/3, —2/3, we have 



|ttGM 7 , b 



1. 7 = —1, the equation l|16fl is (for simplicity we 
choose the changes A — > x, a 7 — - - 
-A/3) 

dx 



dr = 



a/ (a_i + b) x 2 - k 
integrating (|17fl and inverting, we obtain 



(17) 



A(r) 



■ih' 



A - 8vrGM_i 



sinh 



8ttGM_i A 
—3 — ~3 T 



(18) 

The character of this solution is related to the cos- 
mological term A and the curvature parameter k as 
follows, 

(a) For k = 1, the behaviour is inflationary. 

(b) For k = — 1, the behaviour will be inflationary 
if M_i > (A + 3)/8ttG > 0. 

(c) For k = 0, we will solve the original equation 
Ultty. obtaining 



A(r) 



mi exp 



m 2 exp 



2\/H 7r GM_ 1 -|r 



(19) 



here mi and m 2 are integration constants. For 
inflation, the following conditions are neces- 
sary: mi > m 2 and M_i > (A + 3)/87rG > 0. 
This last result generalize that found in [llj . 
and is the same if m 2 = and A = in the 
gauge N = 1. 

2. 7 = —1/3, the equation l|ltj|) is written in the fol- 
lowing form 

dx 



dT = 



with solution 



A(r) 



/8ttGM_i -3k 



|A| 



'1/3 



sinh 



(20) 



(21) 



with A < 0. For inflation, the following conditions 
are necessary M_ 1 > 3k/87tG > 0, implying k = 1 
and |A| > 3M_i. 



3. 7 = —2/3. The equation l|16|l read as 

, dx 

dr = 



^/bx 2 + a_!/ 3 x — k 
with the solution for n — — 1 and A < is 



(22) 



A(r) 



2|A| 



3 3' 



ttGM_ 



^ + -irGM ,) e -^ 
3 3 - M 



—ixGM_ 2 
3 3 



having an inflationary behaviour. 

B. 7 = 1/3, A # 0, any k 

In this subcase, is written as 

AdA 



(23) 



dr = 



|ttGMi - ±AA 4 - kA 2 



with the change of variables u = A 2 , (|24l) is 



r = 



rfi/ 



. /§7rGMl - ±Au 2 - KU 



(24) 



(25) 



whose solutions are, depending on the sign of A 
1. A > 



2AA A 



K 2 + f TTGAMi 



arcsm 



k 2 + fnGAMi 



(26) 



2. A < 



4A Ln 



2W-- | -ttGMi - -AA 4 - kA 2 
V 3 \3 3 3 



-AA - k 



— Ln 



. --CAM - k 



(27) 



C. 7 >0, A = 0, 



For this subcase, Hlfijl is given by 



and introducing the following conformal transforma- 
tion dr = x 37+2 dT, and the change of variable u = 
a-»x~( 37+1 ) — k, the scale factor becomes in the "con- 



formal time" T, 



A(T) 



a 7 (3 7 + l) 2 T 2 



^(3 7 + l)T 



(29) 



which is valid for k < 0. 

When we know the function A(T), we can obtain the 
transformation rule between the times dr and dT, for 
instance 



thus 



dr = A 37+2 dT, 



dr = [^ 7 T 2 - ^ 7 T] 37+1 dT 



(30) 



(31) 



where /i 7 = a "'( 37+1 ^ and v 1 = y/— k(3 7 + 1). 

Now, considering the flat universe (k = 0), we can 
integrate (|31|l . but for consistence between eqs. i|29|l and 
()30|l . we introduce the parameter e in the sense that when 
r = 0, e = T and r = r, e — ► oo 



37 + 2 / , 2(3^ + 1) 

dr = fi 3 T +i / ( x ) 3 T +i ^ x 



After a tedious calculation, we arrive at 



T 



32+2 3(7+1) 

/J. 3t+1 — ; — r - 



3(T+1) 

e 3 7 +i 



37 + 1 
3(7 + 1) 



(32) 



(33) 



3 7 + 1 

Introducing l|33|l into l|29|) we found the scale factor in 
general way 

-3(7 + 1) 



A(r) = 



3 7 + l 



3(7 + 1) 3( 7 + l) 

-r + fj, 2(37 + 1) £ 37 + 1 



(34) 

At this point, we can calculate the behaviour of the scale 
factor for some positive values to the parameter 7 and 
compare them with those found in the standard litera- 
ture, 

1. Dust era, 7 = 0, the scale factor become in the 
gauge N = 1 

2 

A(t) = {[(evrGMo^t + ^V 3 } 3 . (35) 

2. Radiation era, 7 = ~, in the gauge N = 1 
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3. Stiff matter, 7 = 1, in the gauge N = 1 

A(t) = |[247rGMi]H + /i^e-i| 3 . (37) 

Choosing appropriately the parameter e — > 00, we 
obtain the usual results for the scale factor for the 
FRW, with k = and A = Q. 



Also, our method is more general that this employed by 
Faraoni, because when the cosmological constant is in- 
cluded in the formalism, the equation (3.1) in Ref. Q| 
does not reduce to Riccati equation and the procedure 
fails. 



IV. COMMENTS 

The inflationary scenarios for 7 < and matter 
epochs were considered in the FRW cosmological model. 
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